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Abstract 

Using data corresponding to an integrated luminosity of 79 pb" 1 , D0 has 
searched for events containing multiple jets and large missing transverse en- 
ergy in pp collisions at ^/s = 1.8 TeV at the Fermilab Tevatron collider. 
Observing no significant excess beyond what is expected from the standard 
model, we set limits on the masses of squarks and gluinos and on the model 
parameters too and to.j/2, in the framework of the minimal low-energy super- 
gravity models of supersymmetry. For tan/3 = 2 and Aq = 0, with fi < 0, 
we exclude all models with m~ < 250 GeV/c 2 . For models with equal squark 
and gluino masses, we exclude to < 260 GeV/c 2 . 
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Supersymmetry (SUSY) |T| is a symmetry that relates fermions and bosons, and can 
solve the hierarchy problem of the Higgs sector of the standard-model (SM) 0. Minimal 
SUSY extensions of the SM (MSSM) require partners (sparticles) for all standard model 
particles: a scalar partner for each quark and lepton (called squarks and sleptons), and a 
spin-half partner for each of the gauge bosons and Higgs scalars, which form the gluinos and 
the mixed states called charginos and neutralinos. Such models also require the presence of 
two Higgs doublets, and thus four Higgs particles. Each particle in a SUSY model has an 
internal quantum number called .R-parity. If R is conserved, as is assumed in this analysis, 
then sparticle states must be produced in pairs, and each sparticle that decays must contain 
an odd number of sparticles in its decay products. Consequently, in this scenario, the lightest 
SUSY particle (LSP) must be stable, and can thereby provide a candidate for dark matter. 

The most general supersymmetric extension of the SM has over 100 undetermined param- 
eters. Consequently, models have been developed that contain additional symmetries and 
constraints. Some of these involve Grand Unified Theories that include a supersymmetry 
(SUSY-GUTs). In this work, we consider the class of models containing gravity-mediated 
SUSY breaking, called supergravity (SUGRA) models ||. In minimal low-energy super- 
gravity (MLES), the scalar (squark and slepton) masses are unified to a single value m at 
the GUT energy scale, and the gaugino masses are unified to a single value vfi\i%. Three 
other parameters describe the Higgs sector of the model: tan/3, the ratio of the vacuum 
expectation values of the two Higgs doublets; A , a universal trilinear coupling constant; 
and the sign of (x, a mixing parameter in the Higgsino mass matrix. Models in which the 
lightest neutralino (x?) is the LSP, the LSP interacts only weakly. It therefore cannot be 
observed directly, providing an excellent experimental SUSY signature: large missing trans- 
verse energy ($ T ). In such models, squarks (q) and gluinos (g) can decay through a cascade 
of charginos and neutralinos to final states consisting of quarks, leptons, and the LSP. In 
this Letter we describe a search for squarks and gluinos in the jets and lfi T channel. 

The data, corresponding to an integrated luminosity of 79.2 ± 4.2 pb _1 , were collected 
with the D0 detector J| at the Fermilab Tevatron pp collider operating at a center-of-mass 
energy of 1.8 TeV during 1993-95. D0 has three major components: a central tracking 
system, central and forward uranium/liquid-argon calorimeters with towers in pseudora- 
pidity and azimuth of Arj x A<p = 0.1 x 0.1, and a toroidal muon spectrometer. Jets are 
reconstructed using a cone algorithm || with a cone radius of 0.5 in 77 — space. The 
electromagnetic (EM) energy scale is set using the Z — > ee signal. The jet energy scale is 
determined from energy balance in events containing a hadronic jet and a photon candidate. 
$ T is calculated from the vector sum of energy deposited in all calorimeter cells. 

Events were collected using a trigger that required $ T > 40 GeV, and at least one 
calorimeter trigger tower (of size A-q x A0 = 0.2 x 0.2) with transverse energy Et > 5 GeV. 
Offline filtering required $ T > 40 GeV and at least two jets with Et > 8 GeV. 

Backgrounds arise from several sources, among which are it production and the produc- 
tion of W and Z bosons in association with jets. Purely instrumental sources of background 
include QCD multijet production in which jets are mismeasured, resulting in apparent lfi T . 

To remove events with false large $ T , due to detector noise and losses from the ac- 
celerator, we required events to have a summed scalar E T (St) 0.0 < St < 1.8 TeV. 
This requirement has little effect on the efficiency for hard-scattering events with additional 
overlapping soft pp interactions, because these contribute little Et (typical < St > for a 



squark-glunio event is 400 GeV). The position of the primary interaction vertex is required 
to be within 60 cm of the detector center. The initial data sample contains 71,023 events. 
To reject events with large $ T caused by poorly measured jets, we required all jets with 
E T > 15 GeV meet quality criteria based on cluster shape ||, and that the three jets with 
highest Et be within \rj\ < 1.1, or within 1.4 < \rj\ < 3.5. The shape requirements included 
rejecting events in which any jet deposited more than 90% of its energy in the EM portion 
of the calorimeter; this significantly reduced the W — > ev background. 

To select events consistent with being hadronic decays of squarks and gluinos, we required 
at least three jets with E T > 25 GeV, and $ T > 75 GeV. Our trigger was fully efficient 
under these conditions. We also required the leading jet to have E T > 115 GeV, leaving 544 
events. To suppress QCD multijet background we required the azimuthal difference between 
the $ T and a jet of Et > 25 GeV be 5(fi > 0.1, or < {ix — 0.1) radians. To reject events 
where a fluctuation of the second leading jet masks a fluctuation of the leading jet, we also 
required (<50i — 7r) 2 + <50 2 2 > (0.5) 2 , where 1(2) denotes the leading (second-leading) jet in 
Et- To reduce the background from W and Z boson production in association with jets, we 
required Ht > 100 GeV, where Ht is defined as the scalar sum of the transverse energies of 
all but the leading jet. To remove the remaining W —>■ [iv + jets events, we rejected events 
containing isolated muons with p^ > 15 GeV/c. A total of 49 events satisfied all the above 
requirements. 

The average Tevatron luminosity was «9x 10 30 cm" 2 s -1 , and peaked at about 2 x 
10 31 cm -2 s _1 . For the average value, there is a 75% probability of having an additional 
pp interaction accompanying the hard scattering. These additional events contribute many 
charged tracks, which can occasionally cause the soft collision to be chosen as the primary 
interaction vertex, and cause a gross mismeasurement of the $ T . To remove events with 
large ~$ T caused by the misreconstruction of the interaction point, we required the charged 
tracks associated with the central jet of highest Et be consistent with emanating from the 
primary interaction vertex ||. The 15 events passing this criterion formed our penultimate 
event sample. 

The final selection criteria for each (mo, mi/2) point were determined by choosing Ht 
and lp T thresholds that maximized the S/SB ratio, where S is the expected number of SUSY 
events and SB is the combined systematic and statistical error on the background predicted 
from the SM. Table | shows the thresholds used, the number of events expected from SM 
sources, and the number of events observed in the data. 

The largest noninstrumental background arises from the production of it pairs in which 
one t quark decays into jets and the other decays into b£u, where £ = e, /1, or r, and the lepton 
is not detected. We generated it — > btv + jets events using the herwig Monte Carlo [0. 
These were subjected to a detailed detector simulation based on geant [|J and reconstructed 
with the same reconstruction program used for data. We assumed the it production cross 
section of 5.9 ± 1.6 pb |J, which yielded a prediction of 3.1 ± 0.2(stat)ti; 3 (syst) background 
events. 

Comparable backgrounds come from the production of W and Z bosons. Substantial 
~$ T can arise in events with a W boson decaying to leptons where the charged lepton is not 
identified, and in events with Z — > z/z/ or Z — > tt decays. To estimate these backgrounds, we 
generated Monte Carlo samples for W boson events with vecbos |10| (quark hadronization 
simulated using isajet |TTJ), Z bosons with pythia |I2], and WW and WZ events with 



TABLE I. Optimized $ T and Ht thresholds for several regions of MLES parameter space. 
The optimal thresholds were chosen for the specified tuq and m^j values that correspond to the 
listed gluino and squark masses. The next-to-leading order cross sections and the total efficiency 
for signal events, with their combined statistical and systematic uncertainties, the total number of 
events expected from backgrounds, with their statistical and systematic uncertainties, the number 
of observed events, the probability for observing iV bs events or greater given the background 
prediction, and the 95% confidence level upper limit on the cross section are given in the remaining 
columns. Note that the entries in this table are strongly correlated. 
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ISAJET. The detector response was modeled as for the it sample. From all vector boson 
production sources, we predict 2.8 ± 0.8tJx5 events, 85% of which are from W — > Iv and 
Z — > vv decays. 

A major source of background is from events with three or more jets, where one (or 
more) is mismeasured, yielding apparent $ T . To determine this background, we used events 
from 56 pb -1 of data collected with a trigger requiring at least one jet with E? > 85 GeV. 
The trigger was fully efficient for events containing a jet with E? > 115 GeV. Events with 
$ T < 50 GeV were used to determine the instrumental background to events with larger $ T 
using two different estimations. The primary method relied on a Bayesian shape analysis 



We define the quantity D^ = J{8<j>\ — it) 2 + (5(j) 2 — tt) 2 , which has a distribution that is 
strongly peaked at large D^ for events with apparent $ T due to mismeasured jets. For 
it and signal the distribution is less peaked, as shown in Fig. |IJ For the multijet events 
described previously, the shape of this distribution is found to be nearly independent of the 
$ T threshold. To determine the multijet contribution, we performed a three-component (it, 
multijet, and signal) fit to the shape of the D n7T distribution in the data. The backgrounds 
quoted in Table | include the multijet contribution, as determined in this fit. As a check, we 
fit the $ T spectrum of our event sample between 25 and 50 GeV to an exponential in ^ T ; 
extrapolation to higher ~$ T yielded a prediction in agreement with the fit to D^, as shown 
in Table 0. 

To check these background calculations, we relaxed the Ifi T threshold to 50 GeV, and 
obtained predictions of 7.6 ± 0.8Jz|;f events from it and W and Z boson production, and 
35.4 ± 7.9 events from QCD multijet, for a total of 43.0 ± 0.8Jl|;| events from background, 
as shown in Table |. We observed 49 events in the data. 

We note that for all the entries in Table | the number of observed events is greater 
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FIG. 1. Sample D n7T distributions used in the Bayesian shape fitter. Note that the bins 
do not have uniform widths. The top plot shows data passing the analysis requirements with 
]fi T > 75 GeV and H? > 150 GeV (eight candidate events are accepted). The lower plot shows 
the D nn distributions for QCD multijet (thick line), it (dashed line), and an MLES sample with 
mo = 250 GeV/c 2 and mw2 = 60 GeV/c 2 (thin line) for events passing the same requirements. The 
normalizations for the QCD multijet, it and MLES plots are 55.9, 7000, and 350 pb _1 , respectively. 

TABLE II. Comparison of the number of background events expected from QCD multijet 
sources, as obtained from fits to D nn and from extrapolations from lower $ T (see text). The 
uncertainties in the extrapolation do not include the systematic uncertainty due to the dependence 
on the choice of functional form. 
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than the number predicted from background. The results are highly correlated, since most 
rows are subsets of previous rows. The probability for obtaining at least the number of 
events observed for any of the listed cutoffs is more than 2%, and we therefore interpret 
our result as a constraint on the tuq and m\ii parameters of MLES. Simulating squark and 
gluino production and decay with isajet, followed by the same detector response and event 
reconstruction as in our previous simulations, we generated samples at several values of 
m and mi/2, all with the MLES parameters tan/5 = 2, A = 0, and /i < 0. Using the 
next-to-leading order squark and gluino production cross sections from prospino ||14|| , and 



a Bayesian technique with a flat prior for the signal, we determined 95% confidence level 
limits on the parameters. 

Figure [| shows the region excluded by this analysis. Excluded are all MLES models with 



m; 



< 250 GeV/c . For small m , we exclude m~ < 300 GeV/c , and for m 



masses less than 260 GeV/c . In Fig. 
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we show the exclusion contour in the 



m~, we exclude 
(m~, m~) plane 



and compare our results to those from other experiments. We extend significantly the limits 

on squarks and gluinos, especially in the region where m~ > m~. Within MLES models, for 

negative \i and tan/3 = 2, the CERN LEP limits on charginos m~± < 86 to 45 GeV/c 2 pifl 

Xi 

translate roughly to a limit on rrii/2 of 45 GeV/c 2 for small m , and 86 GeV/c 2 for large 
m . Our limit on Tfiyii ranges between 100 GeV/c 2 for small m , and 60 GeV/c 2 for large 
mo- There are no self-consistent MLES models below the solid diagonal line. 
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FIG. 2. The exclusion contour obtained in this analysis (heavy line), the region below which is 
excluded at the 95% confidence level. The thin lines are contours of constant squark or gluino mass 
in the mo - m]/ 2 plane, as indicated. In the shaded region MLES does not contain electroweak 
symmetry breaking, and is excluded apriori. 



In summary, we have searched for events with large Jfi T and multiple jets, and observe no 
statistically significant excess of events beyond expectations from SM processes. This null 
result is interpreted in the context of MLES as an excluded region in the (1710,1711/2) plane. 
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FIG. 3. The limit from this analysis in the (m~, m~) mass plane ( "This experiment" ) . The figure 
also shows curves of previous limits from jets and $ T from D0 (l^Jl^l (hatched) using 7.2 pb _1 of 
data and MSSM parameters tan/3 = 2 and fi = -250 GeV/c 2 , the jets and $ T limit from CDF |l7j 
based on 19 pb _1 of data with tan/3 = 4 and \i = —400 GeV/c 2 (thick dots), the dilepton CDF 
limit [18| from 19 pb _1 of data with MLES parameters tan/3 = 4 and fj, < (dashed-dotted), and 
limits using only direct decays from UA1/UA2 fl|] (dotted) and DELPHI/Mark II [§0| (dotted). 



More recent model dependent CERN LEP limits are given in the text. All limits are at the 
confidence level. The region below the diagonal dashed line is excluded because there the squark 
is the LSP. 
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